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Abstract. Simulation of subgrain growth during recovery is carried out using 
two-dimensional discrete dislocation dynamics on a hexagonal crystal lattice 
having three symmetric slip planes. To account for elevated temperature (i) 
dislocation climb was allowed and (ii) a Langevin type thermal noise was added 
to the force acting on the dislocations. During the simulation, a random ensemble 
of dislocations develop into subgrains and power-law type growth kinetics are 
observed. The growth exponent is found to be independent of the climb mobility, 
but dependent on the temperature introduced by the thermal noise. The in-depth 
statistical analysis of the subgrain structure shows that the coarsening is abnormal, 
i.e. larger cells grow faster than the small ones, while the average misorientation 
between the adjacent subgrains remains nearly constant. During the coarsening 
Holt's relation is found not to be fulfilled, such that the average subgrain size is 
not proportional to the average dislocation spacing. These findings are consistent 
with recent high precision experiments on recovery. 
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1. Introduction 

During the initial stages of annealing of a highly deformed metal, dislocations 
organize into low angle grain boundaries (LAGBs) that enclose cells (subgrains) with 
low dislocation content [lj. At elevated temperatures these walls are mobile leading 
to the gradual coarsening of the subgrains. This phenomenon, known as recovery, 
has a particularly important role, for instance in the nucleation of recrystallization 

m 

The experimental investigation of recovery is rather difficult, since it is obscured 
by the early onset of recrystallization, and usually only a relatively small amount 
of subgrain growth can be observed [21 El [4] . Recently, however, measurements 
were carried out on specimens of certain orientation and deformation modes, where 
the recrystallization is suppressed, and much more subgrain growth occurs [5], El [7]. 
Due to the development of precise orientation measurement techniques [8] , these new 
experimental results shed new light on the physics of recovery. The most important 
findings revealed the abnormality of growth, with large subgrains growing faster 
than small ones |6] , and that the coarsening obeys type 2 kinetics [U [6j [7] . The 
latter means that the time dependence of the average subgrain diameter (D) is 
well-approximated by a power-law 

D n - Dl = ct, (1) 

with .Dq being the cell size at t — and c being an appropriate constant. The growth 
exponent n was found to depend strongly on the annealing temperature and on the 
initial microstructure with values in the range ~ 2 — 7.5 [HE]. Additionally, the size 
distribution of the grains was observed to be close to lognormal, and its variance 
was found to increase faster than its average in time [5j [7] . 

Several methods have been proposed to model subgrain growth. Molecular 
dynamics offers the possibility of investigating the dynamical aspects of recovery 
at the microstructural scale of individual atoms, however, computationally such 
simulations are too expensive [9j HO] . For example in the MD simulations performed 
in Ref. [9] with a volume of 70 x 70 x 1.5 nm 3 and containing 25 subgrains of 15 
nm diameter, the portion of grain/subgrain growth within the timeframe of 7 ns 
is quite limited. Such a simulation is at the limit of both time and length scales 
possible with MD simulations. To overcome these constraints several mesoscopic 
models have been proposed, including Monte Carlo Potts models (e.g. Ref. [TT1[T2] ). 
phase-field method (e.g. Ref. [T3] ) and vertex simulations (e.g. Ref. [HI [T5j US]) - 
all of them relying on input parameters such as grain boundary mobility or energy. 

During recovery mainly LAGBs are moving, where the misorientation between 
adjacent subgrains is less than ^15°. In this case the dislocation cores do not 
overlap, thus it is completely justified to study the properties of LAGBs within the 
framework of dislocation theory (see Ref. [17J or Ref. [18] for reviews). However, 
only a few results have been published concerning dislocation-based modeling of cell 
formation and coarsening. They either focused on LAGB mobility measurements 
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(e.g. Ref. [TH]) or on a limited ensemble (e.g. Ref. |20j). Recently, a new algorithm 
was developed by Bako et al., which, due to a certain multipole method, made it 
possible to study a much larger number of dislocations than before [21, 22j. It was 
shown, that enabling climb leads to cell formation and subsequent subgrain growth 
with the description of this growth given in terms of power-laws [21, 22j. 

In this paper this two-dimensional (2D) discrete dislocation dynamics (DDD) 
technique is developed further. As described in section [2j the modifications are 
undertaken in order to increase the computational precision, include thermal noise 
and use core regularized stress fields to account for dislocation core effects. Due to 
the improvements, the maximum number of dislocations increases to A = 200 000 
in a single simulation run. A method is also developed in this work to obtain the 
corresponding orientation field and the reconstructed subgrain structure. In section 
[3] a detailed analysis of the observed subgrain growth is carried out in terms of 
experimentally measurable quantities. The comparison of the model predictions with 
the experimental findings is given in the second part of section [3] before concluding. 

2. Applied methods 

2.1. Two-dimensional discrete dislocation dynamics 

The simulation of the motion and interaction of individual dislocations is carried 
out using discrete dislocation dynamics in two dimensions. The main features of the 
method are summarized as follows. The system consists of parallel edge dislocations 
which are perpendicular to the plane of the square-shaped simulation area of 
size L x L. Three possible Burgers vector directions ±(cos(m7r/3), sin(ra7r/3)), 
m = 0,1,2 are considered, emulating a hexagonal underlying lattice. Periodic 
boundary conditions are applied and overdamped dynamics are assumed for the 
dislocation motion. 

The stress field of the dislocations is long-range, and therefore at every time step 
all the pair interactions have to be taken into account, which results in an 0(N 2 ) 
algorithm, where A is the total number of dislocations. This time complexity makes 
it impossible to study a large number of dislocations (A > 10 000). To be able to 
consider a larger ensemble of dislocations Bako et al. used coarse grained stress 
fields instead of the analytical ones, which allows the direct calculation of the pair 
interactions to be skipped [21 J. The main idea is to build up a coarse grained discrete 
Burgers- vector density field (a x (l,m),a y (l,m)) defined on an M x M mesh. Let 
ot x {h m ) an d OL y {l, m) denote the x and y coordinate of the net Burgers vector divided 
by the box area in the box indexed by (/,m). From this field the coarse grained 
stress field can be determined on the same mesh using a discrete convolutional 
integral. The approximation here is that the stress field is calculated as if each 
dislocation was at the centre of its cell, so it corresponds to a first order multipole 
approximation. For further details of the method, of the used periodic stress fields 
and of the convolution procedure the reader is referred to [21 J. The most important 
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feature of the method is the reduction of the computational demand which allows 
the inclusion of a much larger number of dislocations than previously (now up to 
200 000 compared to a few thousands). 

In this paper the numerical precision of the 0{N 2 ) algorithm and the more 
advantageous time complexity of the coarse grained method is combined. To this 
end, the following modifications have been included: 

(i) In contrast with the previous method of [21] where, when determining the 
c^(Z,m) fields, each dislocation only contributed to the cell which contained it, 
here the Burgers vector is distributed between the four closest cells weighted 
with a bilinear approximation term. With the notations of figure [TJ if the 
dislocation is positioned between the mesh points (xi,yj) and then 
the weight for the (i, j) node is Wij := (1 — Ax)(l — Ay)/5 2 , for the (i + 1, j) 
node Wi+ij := Ax(l — Ay)/5 2 , etc., where 5 := 1/M. As a result, when a 
dislocation moves from one cell to a neighbour cell, there will be no sudden 
jumps in the c^(7, m) field, and therefore, also in the coarse grained stress field. 

(ii) After performing the convolution, the calculated coarse grained stress field is 
given on the M x M discrete mesh (r cg (i, j)). To obtain a smooth (continuous) 
stress field r™ tp (x,y) we introduce an other bilinear interpolation between the 
cells as 

T ct P ( x >y) '•= ^ijr cg (ij) +w i ^ 1J r cg (i + lJ) +w iJ ^ 1 r cg (iJ + 1) 

+ w i+1J+1 r cg (i + 1, j + 1) (2) 

with the same weighting factors Wij as in the previous point. The same is 
done for the xx and yy component of the stress tensor. The forces acting on 
dislocations are then computed using the continuous stress field r^g tp (x, y). 

(iii) To decrease the inaccuracy introduced by the method a length parameter R e is 
introduced. If two dislocations are closer to each other than i? e , their interaction 
is calculated analytically and not through the coarse grained field. The value 
of R e is chosen to be much larger than the average dislocation spacing. 

(iv) To account for the fact that close to the dislocation core the generated stress 
field does not diverge as 1/r, the core regularized stress fields described in [23] 
are implemented. In this case a new parameter, the core radius r c is introduced 
which is of the order of the Burgers vector. In figure [2] the behaviour of the 
regularized fields is demonstrated by plotting the xy component of the stress 
field induced by a dislocation in its glide plane. 

This modified method is completely continuous, and introduces some 
approximation only at the interaction of far dislocations. To test the level of the 
error, the dependency of the results on R e has been studied. It was found that if 
R e > 2p -0,5 , then this dependence is negligible (here p stands for the total dislocation 
density and p -0 5 is the average dislocation spacing). In the rest of this paper this 
criterion is always fulfilled. 
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Figure 1. The schematics of the bilinear interpolation technique. The dislocation 
position is denoted by the x symbol, and o stands for the grid points of the discrete 
M x M mesh. 



Classical stress field 
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Figure 2. The xy component of the classical and the core regularized stress fields 
in the glide plane of an edge dislocation. According to the double logarithmic plot 
of the inset, the difference between the two curves is negligible after a few r c values 



To mimic the effect of temperature, dislocation climb is allowed using the 
simplest possible mobility rule: 

v g = M g F g and v c = M C F C , (3) 

where v g and v c are the glide and climb velocity of the dislocation, F g and F c are 
the glide and climb components of the acting Peach-Kohler force, M g is the glide 
and M c is the climb mobility. Since the magnitude of the mobilities can be absorbed 
into the time scale [24], the only parameter that affects the dynamics is the ratio 
r] := M c /M g between the two mobilities. 

Besides enhancing dislocation climb, temperature has another important effect 
on dislocation motion: it induces a random movement due to the thermal noise. 
This should be introduced into the model. Usually it is accounted for by adding a 
stochastic component to the force acting on the dislocation segments [25j [26j [27] . In 
this work, when the dynamics are overdamped, it is equivalent to adding a stochastic 
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component to the position (Langevin dynamics). This means that after every time 
step all the dislocations are shifted with Ax g in the glide, and with Ax c in the 
climb directions. The Ax g ^ c values are independent, normally distributed random 
variables with zero mean and, according to the description of Langevin dynamics 
[28] . a half width of 

2 2M g , c k B TAt 

g ' c Al ' [ } 

where T is the temperature, k B is the Boltzmann factor, At is the actual time step 
and Al is the length of the dislocation segment. In this 2D approach it is not 
possible to define A/, therefore, instead of adopting an arbitrary value, in the rest 
of the paper the material parameters are simply merged into a single dimensionless 
effective temperature parameter as 

4vr(l - p)k B T 

TeS •= WKl ' (5) 

where \± is the shear modulus and v is the Poisson's ratio. Note, that the given half 
width of Q is the same as in [27] , but slightly different from the results reported in 



An important note has to be made at this point. If the elastic constants and the 
Burgers vector of for example aluminium at T = 300 K are taken with Al = 10 -6 
m as the typical value for the average dislocation spacing, then T e ff = 1.3 x 10 -4 
is obtained, which is three orders of magnitude lower than the maximum value 
used throughout this paper. However, for real 2D crystals with a triangular lattice 
(the physical system that directly corresponds to the present setup), such as dusty 
plasmas [29j [30] , vortex lattices in type II superconducting films [2U E2], colloidal 
crystals (33], [34] and foams [35], T e $ can be of the order of 1 due to the difference in 
material parameters. So, the difference in the T e ff parameter can be attributed to 
the fact that a 3D system is modelled in 2D. 

The equation of motion ^ is solved by a fifth order Runge-Kutta method [36] . 
Annihilation events are introduced if two or three dislocations with zero net Burgers 
vector are closer to each other than a certain predefined value rf a nnih- 

The presented model is simple and two-dimensional, such that a quantitative 
agreement with the experiments is beyond its capabilities. However, 2D methods are 
rather usual for grain/subgrain growth modelling and provide excellent qualitative 
comparisons with experiments, enabling an insight into the mechanisms at play 
which cannot be deduced from experiment which will be shown in this paper. 

2.2. Simulation parameters and dimensionless variables 

The simulations are started from random distributions of A = 200 000 dislocations. 
The initial number of dislocations is equal on all three slip planes, and the net 
Burgers vector is zero. The core radius is set to r c = L/2000, which is around 
20% of the initial average dislocation spacing, and decreases to around 5% during a 
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typical simulation. For the coarse graining mesh M = 512 is chosen. The parameter 
rj is varied between 0.02 and 0.25, and the effective temperature T e s between and 
0.2. The annihilation distance rf ann ih was set to be equal to the core radius r c . 

According to the equation of motion ([3]), the material parameters can be 
absorbed into the time scale [24]. In contrast to the situation in [24], here the 
dislocation density is not constant and a new length scale r c has been introduced, 
so it is advantageous to use L as the normalizing length scale instead of the average 
dislocation spacing. The dimensionless variables (denoted with prime ( / )) are, 
therefore, defined as 

x' := x/L and t' := tM g Gb 2 /L 2 , (6) 

where b is the magnitude of the Burgers vector and G = /x/[27r(l — v)\. 

2.3. Characterization of the subgrain structure 

During the simulations subgrains form and then grow, in agreement with earlier 
simulations [20], [21] . A typical sequence of dislocation configurations observed (if 
dislocation climb is present) is plotted in figure [31 Although subgrain formation is 




Figure 3. Typical microstructure evolution observed during the simulations. In 
the images t' denotes the normalized simulation time and TV the total number of 
dislocations (which is equal to p' in normalized units). 

apparent, it is not straightforward to define the distinct subgrains from the position 
data of the individual dislocations. The method applied is now described (and is 
similar to the one of |20j). According to the Kroner- Kosevich continuum theory 
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of dislocations, in the present case, the orientation field cj(r) is connected to the 
dislocation density tensor a as [31 



d x uj = a X) and d y uj = a y , (7) 

where the notations a x := a^i and a y := 0^2 are introduced{|] These quantities are 
the continuum versions of a x ^ y (l,m) defined in section |2T] Equations ([7| can be 
transformed into a Poisson equation 

Auj = d x a x + dyoty, (8) 

which can be solved efficiently on a discrete grid iteratively by a multigrid method 
[36] . Consequently, periodic boundary conditions are automatically fulfilled. Thus 
from the dislocation positions the discrete a x ^ y (l,m) fields can be constructed, and 
then u)(l,m) follows from This procedure can be performed on an arbitrarily 
smooth K x K mesh. Figure [4](b) shows the orientation map corresponding to the 
configuration of figure |i](a) for K = 512. It is important to note that without 
the introduction of core regularized stress fields it is not possible to define b in 
a straightforward manner, and therefore, the 6 ~ b/d misorientation cannot be 
defined. In the present setup, however, one can assume b := r c , leading to defined 
uj values in figure [3](b). So, the misorientation values are directly comparable to the 
experiments. 

The subgrain structure is then derived by locating close to constant orientation 
regions on the uj field. Figure [4](c) shows the structure derived from figure [3](b), 
the colours referring to the average orientation in the subgrain. According to 
figure [4](c), the large subgrains are determined confidently, while there is some 
arbitrariness in the definition of the small ones. This is the direct reflection of 
the unavoidable uncertainty of the subgrain definition. After the subgrain structure 
has been reconstructed the average area A and the standard deviation 5 A can be 
measured. The linear size is then D := \T~A and 8D := ^/6A. 

The subgrain sizes are of course not uniform, they obey a certain distribution 
P(A), where A denotes the area of a subgrain. From the reconstructed structure of 
figure [4](c) P(A) is easily obtained. For better numerical accuracy the cumulative 
distribution function F(A) of P(A) was constructed and then averaging was 
performed over statistically equivalent realizations of the system. In figure [4](d) 
F(A) is plotted, and one can conclude from the fit that the distribution is close to 
lognormal (just like in experiments [5j Ej). Note that here and in the rest of the 
paper, ensemble averaging is always performed over 16 independent simulation runs. 

J In the Kroner- Kosevich continuum theory a quantity called the 'dislocation part of the relative 
rotation' is introduced, which is the non-elastic component of the rotation field, and is denoted by 
0. The introduced uj field corresponds to the z component of this field. 
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Figure 4. The identification of the subgrain structure, (a) A sample dislocation 
configuration, (b) the corresponding orientation map (in degrees), (c) the identified 
subgrains (the colours refer to the average orientation within the subgrains) and 
(d) the averaged subgrain size distribution. The continuous curve is the measured 
distribution and the dashed one is the fitted lognormal distribution. 



3. Results 

3.1. Power-law growth 

First the effect of the climb mobility is investigated, whereby simulations with 
different rj values have been performed at T e $ = 0. As seen in figure [5| the average 
subgrain size D f follows power-law kinetics. Since the smallest identifiable subgrain 
is limited by the resolution of the mesh on which ^ is solved, at small i! values a 
higher average is obtained than the real value. On the other hand, when D' ~ 0.1, 
the size of the largest subgrains is already in the order of the system size L, which 
biases the average in the other direction. Nevertheless, the power-law regime lasts 
for more than 1.5 decades, which confirms that the coarsening obeys type 2 kinetics. 
From the fit n = 2.9 ± 0.2 is obtained, which is not far from the experimentally 
observed low-temperature values [TJ [6] . 

The other consequence of figure [5] is that the power-law exponent does not 
depend on rj. More precisely, a slight deviation is only seen when r] is larger than 
0.1. This is consistent with the results of Hartmaier et al. [38] who studied the 
creep properties of thin films. Therefore 77 = 0.1 is implemented in the rest of the 
work presented in this paper. Moreover, the curves almost overlap when they are 
plotted versus t'rj indicating that the climb rate only modifies the time scale of the 
simulation but not the nature of the dynamics. 

The constancy of the power-law exponent contradicts experimental findings 
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Figure 5. The average subgrain size evolution at T e ff 
mobilities. The values were plotted against t'rj. 



and different climb 



which report decreasing n with increasing temperature [HE]. However, the addition 
of thermal noise to the dislocation motion (as described in section 2.1) leads to the 
expected increase of the power-law exponent, that is, to the decrease of n. This is 
demonstrated in figure [6j where it is seen that the exponent of n = 2.9±0.2 decreases 
to about 2.2 ± 0.2. The numerical values of the exponents are different from the 
experimental values, which is one of the limitations faced by a 2D model studying a 
3D phenomenon. In addition, several mechanisms that hinder boundary motion, 
such as dislocation-point defect interactions, are not considered. The tendency 
observed, however, is in good agreement with the experiments [HE]. 




Figure 6. The average subgrain size as a function of time at 77 = 0.1 and different 
T e ff values. 



The exponent was found to depend on the initial number of dislocations N as 
well as on the core radius r c . As seen in figure [7| n increases as the core radius 
and the initial density grow and as expected, only the r c /po°' 5 ratio, i.e. the ratio 
between the core radius and the average dislocation spacing, is important. A possible 
reason for this dependence could be that the mobility of a dislocation wall increases 
with decreasing core radius. However, for a thorough explanation further studies 
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Figure 7. The average subgrain size as a function of time at r] = 0.1 and different 
r c and p f values. 



It is also instructive to look at the evolution of the average misorientation 
between adjacent subgrains. Starting from the reconstructed subgrains of figure [i](c) 
the mean uncorrelated misorientation 9 U was measured which is by definition the 
mean misorientation between all possible neighbour subgrain pairs. As seen in 
figure [8| a slowly growing average misorientation was found which is non-sensitive 
to the climb mobility and its growth rate decreases with increasing temperature. In 
experiments 9 U was found to be nearly constant [6j . In this work the growth rate is 
low, and is even decreasing at high temperatures, thus providing good qualitative 
agreement. 
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Figure 8. The average misorientation as a function of the average subgrain size. 



The discussion now returns to growth exponents. It was found above that 
the exponents n are slightly below the experimentally obtained values. They are 
surprisingly high, however, if compared to the predictions of the existing models 
of subgrain growth kinetics. Models based on dislocation climb predict parabolic 
subgrain growth, i.e. n = 2 [H SI 1321 40j. As a possible solution for the discrepancy 
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it was suggested that the thermally activated migration of the ledges is the rate 
controlling mechanism [4], or the gradual decrease of the average misorientation is 
responsible for the increased exponent [6] and not the climb. According to these 
simulation results it is evident that a simple climb model is able to account for an 
exponent of n = 3.3, even with a slightly increasing average misorientation. Thus it 
is likely that some important feature of subgrain growth is missing from the existing 
climb-based models, which should, therefore, be revised. 



3.2. Discontinuous subgrain growth 

During discontinuous (or abnormal) subgrain growth large cells grow faster than 
small ones, which leads to an inhomogeneous subgrain distribution during annealing. 
This type of growth was confirmed by recent experiments [5j [6|. A simple 
consequence of this behaviour is that the subgrain structure is not self-similar, and 
the scaled size distribution of subgrains broadens with time. 

From the reconstructed subgrain structures the average linear size of the 
subgrains (D ) and its standard deviation (5D') were measured as described in 
section 2.3, In figure [9] the evolution of the relative scattering 8D f /D' is plotted 
at various simulation parameters. It is evident that the growth is discontinuous, 
and that neither the temperature nor the climb rate has any effect on it. 
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Figure 9. The relative standard deviation of subgrain sizes versus average 
subgrain size. 



3.3. Length scales in the micro structure 

It is commonly assumed that in large dislocation systems there is only one length 
parameter: the average dislocation spacing p -05 . This means that dislocation 
structures are scalable, which in the case of cell sizes is expressed by Holt's relation 

mi 

D = Kp~ - 5 , (9) 
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where K is a suitable constant. Relation ^ was proven experimentally for example 
on deformed iron [41] . on crept LiF [42] or more recently on doped GaAs [43J. It 
is surprising, therefore, that Holt's relation ^ contradicts the observed abnormal 
nature of subgrain growth. Namely, if there is only one length parameter (the system 
is scalable) only self-similar, i.e. normal growth, is possible. But in the present case 
to describe the size distribution at least two length parameters are needed (its mean 
and its deviation) since the observed relative scattering of cell sizes increases, and 
therefore the growth rate cannot be normal. 

It was already mentioned that in these simulations together with the core 
regularized stress fields a new length parameter, the core radius r c was introduced. 
Since Holt's scaling argument only assumes a 1/r type dislocation stress field, this 
could solve the contradiction. Hence, the simulations were repeated with different 
core radii, and the average cell diameter D dependence on the average dislocation 
distance p -0 5 was investigated. Figure 10 shows that the effect of r c is negligible on 
the D — p -0 - 5 relation. This is understood as the average distance of dislocations in 
the LAGBs was found to be at least 3 times larger than the core radius r c , even at 
the highest value of r c . So, it is not the introduction of core regularized fields that 
breaks Holt's relation. 
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Figure 10. The average cell size (D') versus the average dislocation spacing 
(p -0,5 ). Holt's relation would imply a linear dependence. 



The contradiction can be removed if the initial dislocation density rather than 
the instantaneous dislocation density is considered in the scaling relation ([9]). In 



figure [11] the average misorientation 9 U was plotted as a function of the average 
subgrain size D at different initial dislocation densities p - The latter simply 
indicates the different number of initial dislocations in the L x L simulation box. In 
addition, the misorientation was divided by the initial average dislocation spacing. 
The curves fully overlap revealing 9 U = Kp^ 0-5 , with K depending on time, average 
subgrain size, etc. 

It is reasonable to assume that during the initial stages of the simulations, 
when the first dislocation walls are formed from the random arrangement of 
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Figure 11. The average misorientation weighted by the initial average dislocation 
spacing as a function of the subgrain size. The lower is the initial density, the lower 
is the observed misorientation. 



dislocations, the average distance of the dislocations within the walls must scale 
with Pq ' 5 . According to this discussion, this distance, which is directly related 
to 9 U through Frank's formula, remains nearly constant throughout the growth 
process. The slight increase observed can be related to the fact that the mobility 
of LAGBs is proportional to the misorientation [I], therefore, LAGBs with higher 
misorientation wipe out the sparse ones (having lower misorientation). The picture 
which emerges is in complete agreement with the experiments, where the constancy 
of the misorientation is valid even for specimens with quite different initial 9 U [6|. 

Note that contrary to the present investigations, in the previous work of Bako 
et al. [21], Holt's relation was found to be satisfied throughout the growth process. 
After the more precise analysis of this paper, the conclusion obtained is that those 
results were due to numerical noise. The latter is similar to thermal noise, such 
that those results should be regarded as if very high temperature was applied (see 
figure [6]) . 

The next question to be addressed is why Holt's relation is valid in a wide 
range of systems with cell structure. The possible reason is that in those cases 
subgrain formation is induced by external stress, rather than being a simple 
growth phenomenon. It is well-known, that during creep, the subgrain structure is 
characterized by a dynamic equilibrium, that is, the plastic strain is produced by the 
moving subgrain walls, but still, the average subgrain size remains constant [44J. In 
this state, the net dislocation annihilation must be balanced by creation mechanisms. 
The latter is absent in this model and also during recovery. In summation, from 
this investigation it is evident that dislocation creation is the key to the dynamic 
equilibrium during creep, and to the fulfilment of Holt's relation. 
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4. Conclusions 

Cellular dislocation patterning and its growth phenomenon has been investigated 
within a simple 2D discrete dislocation dynamics model. The effects of the climb 
rate, the thermal noise and the size of the dislocation core radius on the kinetics of 
growth have been studied. The main results of the paper are as follows: 

(i) The coarsening follows power-law kinetics, with a growth exponent not 
depending on the dislocation climb mobility but rather on the thermal noise 
and core radius. 

(ii) The size distribution of the subgrains was close to lognormal, which broadened 
with time. This means that growth is abnormal (discontinuous). 

(iii) Holt's relation is not fulfilled in the conventional way - the average subgrain size 
is proportional to the initial dislocation spacing not to the temporal spacing. 

(iv) The average misorientation between the adjacent subgrains slightly increases 
throughout the growth process. 

(v) Despite the simplicity of this model, all of the above results are in qualitative 
agreement with the experiments. 

(vi) In contrast to previous dislocation climb-based models, our simulations yield 
growth exponents considerably larger than n = 2. 

These results indicate the success of using 2D discrete dislocation dynamics in 
modelling subgrain coarsening. 
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